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ABSTRACT 

We analyzed the relation between the corotation radii and the galactic radii at 
which breaks or changes of slope of the metallicity gradients occur in spiral galax- 
ies. With this purpose we compiled the results from the literature on rotation curves, 
corotation radii and radial metallicity distributions of 27 galaxies, of which 16 were 
considered qualified to be studied in the context of this work. We re-scaled all ref- 
erences of each galaxy to a same framework in order to compare the results and to 
identify the radii where breaks and changes of slopes are found, when non-linear mod- 
els fit the radial metallicities better than a linear model. In most galaxies we have 
found minima and breaks in radial metallicity near the corotation radius, revealing a 
significant correlation between these two radii, as it occurs in our Galaxy. The results 
are interpreted as a consequence of long-lived spiral structures, in which the star- 
formation rate depends on the distance to the corotation radius, producing secular 
effects in the observed radial metallicity distributions. 
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1 INTRODUCTION 

The mechanisms that regulate the star-formation rate and 
the building up of the chemical abundances in disks of 
galaxies, and in particular, the role of the spiral structure, 
are not yet fully understood. In the past, most stud- 
ies of radial metallicity distributions, especially of the 
oxygen abundance distribution, re vealed monotonically de- 
creas ing 
19921: 



r adial gradients (eg. Vila-Costas fc Edmunds! 



iKennicutt. Oev. Zaritskv fc Huchraf 



Zaritskv, Kennicutt Jr. fc Huchral 1994 ) 



19931 ; 
This simple 

behaviour of metallicity gradients was easily repro- 
duced by chemical evolution models in which the 
most important parameters are supposed to 



most important parameters are s upposed to vary 
smoothly with radius (see eg. Molla. Ferrini fe Diadll996l ; 
IChiappini. Matteucci fc Grattor] 1997). For these reasons 
the linear behaviour of the gradients (in logarithmic scale) 
seemed to be universal. However, the presence of breaks in 
the oxygen abundance slopes had already been observed in 
a number of galaxi es in the decade of 1990, for instance in 
M101 (NGC5457) (Zaritskv 1992; I Vila-Costas fe Edm und's] 
1 19921 ). NGC5194 ijVila-Costas fc Edmunds! Il992i ). NGC 
3319 and M83 (NGC 5236) (|Zaritskvl 1 19921 ). NGC3359 
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jMartin fc Rovl Il995l ). NGC 1365 l|Rov fc Walshl Il997l ). 
The rea l exist ence of these breaks was contended by 
Pilyugin (2003), who claimed they are due to a systematic 
error involving the excitation parameter, which would 
overestimate the oxygen abundance in low excitation H H 
regio ns. But at least i n a few cases, like M83 (jBresolin et al.l 



120091). NGC3359 by IZahid fe Bresolm] l|2011l) . NGC 4625 
( Goddard et~alll201ll ), new studies or re-examination of the 
gradients with new metallicity diagnostics showed that the 
breaks indeed exist. 

In general the breaks in abundance gradients are at- 
tributed to variations of gas density, the presence of a bar, 
or coi ncidence with the corotation radius (Rcr)- In a recent 
work, IScarano Jr.. Lepine fc Marcon-Uchidal (1201 ll ) investi- 
gated 3 galaxies for which the Rcr was supposed to be in 
the visible part of the spiral structure. They found evidences 
for variations of gradient slopes, and argued that they are 
associated with the corotation resonance. 



In the case of our Galaxy a long debate con- 
cerning the existence of breaks in the gradient also 
occurred. After many works report i ng lin ear gradients 
iTwarog. Ashman fc Anthonv-Twarod (|l997l ). based on a 
sample of Open Clusters, were the first to show that there 
is an abrupt step down of 0.3 dex in the [Fe] abundance, at 
about 1 kpc beyond the orbit of the Sun. However, several 
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authors, analyzing similar sample s of clusters, did not con- 
firm th e existence of the s t ep (e g. Rolleston et aljfeooch . In 
parallel lAndrievskv et al performed abundance 

measurements of several elements in Cepheids, and found 
that the metallicity gradient presents a steep slope around 
5-6 kpc from the center, a plateau from about 6.5 kpc to 
10 kpc, followed by a new decreasing gradient at about 10 
kpc (quoting here the radii given in the original papers). 
iLepine et all (|201ll ). reanalyzing the most recent data on 
Open Clusters and Cepheids, confirmed the existence of 
a sharp step in the abundance gradient, and showed that 
it coincides with Rcr, which is now accurately known in 
our Galaxy. The theoretical explanation proposed for the 
metallicity step in our galaxy is the "pumping out" ef- 
fect of corotation, which produces gas flows in opposite 
directions on the two sides of the resonance. This origi- 
nates the ring-shaped gap in the distributio n of gas ob- 
served by lAmores. Lepine fc Mishurovl (|2009ri that isolates 
the inner and outer sides of the gap one from the other. 
As a consequence, the metallicity has evolved in an inde- 
pendent way on the two sides. A rec ent study of NGC1365 
l|Elmegreen. Galliano fc Alloinl 120091 ') observed the flow of 
gas moving out from corotation, in opposite directions on 
each side. 

The existence of plateaux, sharp steps, or rising metal- 
licity gradients in the outer regions of galaxies, should not 
be regarded as minor effects that are superimposed on a 
general, well understood, linear behaviour of metallicity gra- 
dients. These features reveal that important processes like 
radial gas flows do exist and should be taken into account 
in chemical evolution models. 

In order to check if the same effect observed at the coro- 
tation resonance in our Galaxy plays a major role in deter- 
mining the shape of the metallicity gradients in other spiral 
galaxies, we analyzed the data from more than a hundred 
references on metallicity distributions, rotation curves and 
corotation radii (or their equivalent spiral pattern speeds) 
to verify the existence of a correlation between the radius 
of corotation and the radius where changes in slopes, or 
even of discontinuities, in the metallicity gradients occurs. 
In Section 2 we examine the case of our Galaxy, with the 
purpose to establish the strategy to locate minima, points of 
inflexions and/or discontinuities. In Section 3 we introduce 
the data compiled from the literature and the procedures to 
standardize the results from different references. In Section 
4, we explain the method used to locate the breaks, minima 
or inflexions in the metallicity distributions. In section 5 we 
present the correlation between the metallicity breaks radii 
and corotation radii. The remarks on individual galaxies are 
collected in Section 6. We conduct the discussions in Section 
7, and finally, in Section 8 we summarize the main results 
and discuss some of the consequences of the present work. 



2 THE ABUNDANCE GRADIENT 

DISCONTINUITY IN OUR GALAXY 

It is important to fully understand the connection between 
metallicity features and corotation in our Galaxy before 
proceeding to the analysis of other galaxies. The coro- 
tation radius of our Galaxy is accurately known. Many 
papers in the past pointed that corotation was situated 
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Figure 1. Two examples of abundance gradients in the Galac- 
tic disk, from the literature: (a) The [Fe/H] distribution of Open 
Clusters (Lepine et al. 2011), (b) The [O] abundance of OB stars 
from Daflon et al.(2004). In Panel (b) the galactoccntric radii 
of the sources where re-calculated for Rq= 7.5 kpc. The oxygen 
abundance scale was normalized to the solar abundance using 
[0/H]q = -8.6. In the two panels the dashed lines are fits of two 
lines (see text), the vertical line indicates the position of corota- 
tion, and the thin solid curve the fit of the data by a polynomial 
with the degree that best fits the observations. 
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others). More recently a few direct methods not only con- 
firmed these previous estimates but also allowed to decrease 
the uncertainty on the results to a few hundred parsecs. 
By 'direct' we mean that we exclude methods like chemical 
evolution models, n-particle simulations, etc, the results of 
which strongly depend on a number of hypotheses. Among 
the direct methods we admit the determination of the pat- 
tern speed by integrating the orbits of open clusters back 
to their birthplaces. The result, R C r = (1.06 ± 0.08)7? Q , 
does not ha ve any strong depende nce on the rotation curve, 
as shown bv lDias fc Lepina (|2005T ). Note that a contestation 
of this result, claiming that different arms have different ve- 
locities, was refuted bv lLepine et al.l l|201ll ) (L+8 hereafter) . 
Other direct observati ons are the position of the ring-shaped 
gas void at corotation (jAmores et al.ll2009h , and the position 
of the square-shaped spiral arm associated with the 4:1 res- 
onance (L+8). These observations of resonances point to a 
slightly larger Rcr, but still within the error bar of the de- 
termination of Rcr mentioned above. We adopt here Rq = 
7.5 kpc and Rcr= 8.1 kpc. The choice of the short scale 
of the Galaxy instead of -R© = 8.5 kpc recommended by the 
IAU is based on many recent measurements, as discussed 
in detail by L+8. However, the conclusions reached in this 
paper do not depend on that choice. 

We illustrate in Figure [1] the behaviour of two different 
tracers of chemical abundance in the region of corotation of 
our Galaxy. In Panel [TJl we show the Fe abundance of open 
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clusters as a function of radius, taken from L+8. There is 
a clear step down in metallicity and a gap in the distribu- 
tio n of objects at 8.5 kpc. This step was first discovered 
bv lTwarog et all (|l997T ). but its connection with corotation 
was not noticed at that time. The distributions of metallic- 
ity shown in the figure are satisfactorily described by two 
plateaux, indicated by the vertical dashed lines. There is 
possibly some overlap in galactic radius between the high 
level and low level metallicity plateaux. Whereas this is one 
of the clearest evidences for the existence of a metallicity 
discontinuity, it is not the best choice for a comparison with 
other galaxies. In extragalactic studies the focus is on O 
abundance in H II regions, while the Fe gradients are not 
available. Furthermore, Open Clusters are objects older than 
H II regions, and their age give them the possibility of mi- 
grating out from their birthplace radius, as discussed by 
L+8. For instance the clusters situated beyond 13 kpc are 
old clusters that were probably born at a smaller galactic 
radius. Therefore, the extension of the low-level plateau to 
large radii probably does not correspond to a real metallicity 
plateau of the local gas. Similarly, the large radial overlap 
between the high-level and low-level plateaux could be ex- 
plained by the radial migration of the clusters. 

In Panel [Ha we show the ox ygen abundance data of OB 
stars of lDaflon fe Cunhal (|2004f ). This is the best example of 
Galactic abundance gradient to be compared with those of 
other galaxies, since the objects are quite young (less than a 
few Myrs, which excludes the possibility of migration) and 
oxygen is the element being measured. Two straight lines fit 
well the data, and the step of the gradient is quite clear. Be- 
yond the corotation, the gradient presents a positive slope. 
As it will be shown in later sections, such positive slopes are 
not rare in other galaxies as well. 

In the two panels of Figure [Tj we also fitted the data 
with polynomials of different degrees (we selected the de- 
gree which gives the best fit, according to the discussion in 
Section 4). It is remarkable that the fitted polynomials show 
a point of inflexion that coincides with Rcr- This analysis 
suggests that when a break or step of the gradient is clearly 
visible, the best approach is to adjust two independent lines, 
one on each side of this transition point. However, when the 
data presents only smooth variations of the slope, the min- 
ima of a fitted polynomial can be a good indicator of the 
position of corotation. 

We should mention that there are in the literature on 
our Galaxy many works presenting the abundance of some 
elements as a function of radius. The results are based on 
all sorts of methods and objects. Only rarely one can clearly 
distinguish the discontinuities in the metallicity at Rcr- The 
main reasons are either too few data points, large error bars 
in the distances or in the metallicities, or the use of old 
objects which had time to migrate. 



[T] contains the basic data for these galaxies; the references 
are given in Tables [ATI and [A"2l in the Appendix. 

As expected, the measurements made by many groups 
of research and at different epochs are based on a variety 
of methods and assumptions. For this reason the following 
criteria were used to select the data and to standardize the 
measurements: 

Geometrical Parameters: The distances of objects to the cen- 
tre of a galaxy depend fundamentally on the distance to the 
galaxy and on projections parameters of its disk on the plane 
of the sky (inclination and position angle). The first step 
of this work was to recalculate the radial distances of the 
objects given in different papers, adopting for each galaxy 
the distance and projection angles that we determined. We 
compared the results obtained by the traditional method of 
fittin g ellipses to isophotal curves, and from spiral arms fit- 
ting (|Mall200ll ) in images available in DSS. For the present 
sample of galaxies the results are in agreement with each 
other and also with the results on the literature, especially 
those derived from velocity fields (see references in the ap- 
pendix). Since the differences of the resulting angles based 
on the different methods are smaller than 15°, we adopted 
the median of all results to derive the values in Table [T] 
The most critical discrepancies (of almost a factor 2 in some 
cases) were on the distances to the galaxies. We selected the 
distances given in the literature according the method used 
to determine them. The preferences among the methods fol- 
lows the sequence: 1-) Cepheids; 2-) Supernovae; 3-) Hubble 
Flow, adopting H=73 km/s/Mpc. 

Metallicity Distributions: The most common way to measure 
the metallicity distribution in spiral galaxies makes use of 
the nebular lines from H II regions. In a strict sense, "metal- 
licity" refers to the abundance of elements heavier than He- 
lium. Among these elements oxygen is the most abundant 
and easiest to be measured. It is also the element which 
has abundances available for the largest number of galaxies. 
Consequently, we concentrated our study on this element. 
Since O is mainly released by Type II Supernovae, objects 
that are strongly associated with the spiral arms of galax- 
ies, it is also a convenient element to investigate the effects 
of the corotation resonance on the radial metallicity distri- 
butions. There are several calibrators that allows calculate 
the oxygen abundances based on the ratio of line fluxes, 
going from the most direct methods, based on temperature- 
sensitive line ratios (Te- method hereafter), to some sort of 
statistical methods that do no use oxygen line fluxes. Since 
each calib rator may be af fect ed by different sys t emati c ef- 
fects (see iPilvugini <|2003h or iKewlev fe Ellison! (|2008h for 
example). For this reason we propose to study the changes 
in the radial metallicity distributions of a galaxy using all 
calibrators available and verifying the consistence between 
the results. 



3 DATA SAMPLE AND NORMALIZATION 

The data used in this work were compiled from the litera- 
ture. We selected galaxies for which there were one or more 
papers giving their corotation radii, radial metallicity distri- 
butions and rotation curves. This criterion was fulfilled by 27 
galaxies with their data distributed in 130 references. Table 



Rotation Curves: It is not uncommon to find differences on 
the rotation curve of a galaxy in different papers, especially 
when they were observed with different techniques. Prefer- 
ence was given to: 1-) Rotation curves composed by radio 
and optical observations; 2-) when optical and radio observa- 
tions are available, but they have been presented separately, 
we joined them to use the resulting curve; 3-) rotation curves 
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Table 1. Selected galaxies and their projection parameters (see details in the text). The symbols used are i for inclination and PA for 
the position angle, both in degrees, D is the distance to the objects, in Mpc, and s is the scale (kpc/arcmin) for the deprojected galaxy. 
The indices between brackets are the references for the distance (D) and the rotation curve (V(r)) used along the paper. 



Galaxy 


AR 


Dec 


i 


PA 


D 


s 


Method 


References 




J2000 


J2000 


[°] 


(J2000) [°] 


[Mpc] 


[kpc/'] 




D 






IC0342 


03:46:48.5 


+68:05:46 


21 


130 


3.3 


0.96 


Cph 


[100] 




[113] 


NGC0224 (M31) 


00:42:44.3 


+41:16:09 


74 


36 


0.73 


0.21 


Cph 


[126] 




[113] 


NGC0598 (M33) 


01:33:50.9 


+30:39:36 


53 


25 


0.8 


0.23 


Cph 


[67] 




[113] 


NGC0628 (M74) 


01:36:41.8 


+15:47:00 


16 


25 


6.7 


1.95 


SNII 


[127] 


[23] 


+ [58] 


NGC1068 (M77) 


02:42:40.7 


-00:00:48 


32 


72 


15.7 


4.57 


Hub 


[24] 




[113] 


NGC1232 


03:09:45.5 


-20:34:46 


20 


269 


21 


6.11 


Hub 


[24] 




[120] 


NGC1365 


03:33:36.4 


-36:08:25 


46 


222 


18.3 


5.32 


Cph 


[111] 




[113] 


NGC1530 


04:23:27.1 


+75:17:44 


50 


30 


35.8 


10.41 


Hub 


[24] 




[72] 


NGC1566 


04:20:00.4 


-54:56:16 


16 


60 


18.2 


5.29 


Hub 


[24] 




[86] 


NGC2403 


07:36:51.4 


+65:36:09 


60 


123 


3.47 


1.01 


SNII 


[128] 




[113] 


NGC2543 


08:12:57.9 


+36:15:17 


60 


46 


33.5 


9.74 


Hub 


[24] 




[72] 


NGC2903 


09:32:10.1 


+21:30:03 


69 


21 


6.52 


1.9 


Hub 


[9/11 
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NGC3031 (M81) 


09:55:33.2 


+69:03:55 


60 


157 


3.63 


1.06 


Cph 


[41] 




[113] 


NGC3319 


10:39:09.4 


+41:41:12 


58 


42 


14.3 


4.16 


Cph 


[101] 


[79] 


+ [48] 


NGC3359 


10:46:36.8 


+63:13:27 


51 


172 


14.4 


4.19 


Hub 


[134] 




[3] 


NGC3621 


11:18:16.5 


-32:48:51 


65 


345 


6.6 


1.92 


Cph 


[42] 




[7] 


NGC4254 (M99) 


12:18:49.6 


+14:24:59 


11 


68 


32.3 


9.4 


Hub 


[2d] 




[64] 


NGC4321 (M100) 


12:22:54.9 


+15:49:21 


27 


146 


16.1 


4.68 


Cph 


[39] 




[113] 


NGC4736 (M94) 


12:50:53.0 


+41:07:14 


50 


125 


4.93 


1.43 


Hub 


[24] 




[113] 


NGC4826 


12:56:43.7 


+21:40:58 


59 


114 


5.48 


1.59 


Hub 


[24] 




[94] 


NGC5033 


13:13:27.5 


+36:35:38 


63 


170 


12.7 


3.69 


Hub 


[24] 




[113] 


NGC5055 (M63) 


13:15:49.3 


+42:01:45 


57 


107 


15.7 


4.57 


SNIa 


[107] 




[113] 


NGC5194 (M51) 


13:29:52.7 


+47:11:43 


20 


170 


8.9 


2.59 


SNII 


[26] 




[113] 


NGC5236 (M83) 


13:37:00.9 


-29:51:56 


24 


45 


1.5 


1.31 


Cph 


[117] 




[113] 


NGC5457 (M101) 


14:03:12.6 


+54:20:57 


18 


39 


7.4 


2.15 


Cph 


[57] 




[113] 


NGC6946 


20:34:52.3 


+60:09:14 


44 


58 


3.76 


1.09 


Hub 


[24] 




[113] 


NGC7479 


23:04:56.6 


+12:19:22 


46 


41 


34.9 


10.15 


Hub 


[24] 




[49] 



from radio velocity fields only; 4-) rotation curves from op- 
tical velocity fields only; 5-) rotation curves from long slit 
data. All rotation curves were recalculated considering the 
distances, inclinations and position angles we adopt in this 
work (Tabled]). 

Corotation radius and/ or spiral pattern speed: These two pa- 
rameters are easily converted from one to the other if we 
know the rotation curve (cj p = V(Rcr)/Rcr)- In principle 
we did not make any distinction between methods to eval- 
uate the corotation radii or the spiral pattern speeds. How- 
ever, we confront the values obtained in the literature with 
the constraints imposed by the rotation curve and the ex- 
tension of the spiral structure. When several measurements 
where available in the literature, we only considered the me- 
dian of the values which are not in conflict with these con- 
straints. In our effort to standardize the radii of corotation 
we faced two situations: 1-) cases when there are differences 
between our standard distance to the galaxy and the dis- 
tances adopted by the consulted references, but similar ro- 
tation curves were considered; 2-) cases where there is a sig- 
nificant difference between our standard rotation curve and 
the rotation curve from the references on corotation. In the 
first case we just scaled the corotation radius or the spiral 
pattern speed to the new distance. In the second case it is 
necessary to convert the pattern speed of the original refer- 
ence to corotation radius, using the original rotation curve 
and then scale this radius to the adopted distance. The final 



corrected pattern speed is obtained by converting back the 
scaled corotation radius using the standard rotation curve. 
The uncertainty limits (or extremities of the error bars) of 
the original paper are converted to the new system using the 
same procedure. Since rotation curves are usually not lin- 
ear, this conversion produces some asymmetric error bars. 
In these cases we adopted the largest values of uncertainties. 



4 METHODS TO IDENTIFY BREAKS IN THE 
METALLICITY DISTRIBUTION 

4.1 Preliminary tests with Monte-Carlo 
simulations 

We have seen from the example of our Galaxy, that in some 
cases, when there are plateaux separated by a step, one be- 
fore and one after the corotation radius, an inflexion point 
reveals the radius of the break. However, we also found cases 
of galaxies for which a change in the gradient slope is clearly 
associated with corotation, without a noticeable step in the 
level of metallicity. In such cases a good representation of 
the data can be obtained by fitting them with two lines 
connected at their intersection. There are still other cases 
where there is a minimum in the metallicity at the corota- 
tion radius, because the slope of the gradient changes from 
negative to positive. In such cases a polynomial fitting is a 
convenient procedure to identify the minimum. 

A question that arises is: how can we check that the 
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function that we adopt to represent the abundance distri- 
bution is a better choice than a straight line, which is the 
usually accepted fit? 

In order to establis h the significance o f a m odel with 
a break in the gradient, IZahid fc Bresolin] (|201ll ) has per- 
formed likelihood ratio F-test to the galaxies they studied, 
but we did not adopt the same approach. In many cases it 
is not possible to standardize the data in the literature in 
reason of the lack of information on the flux of some strate- 
gic lines. Another problem arises mainly from the different 
ways the uncertainties are evaluated in the literature, and 
the F-test is particularly sensible to these differences. This 
problem would force us to discard many objects in a sample 
that it is already small, so we opted to follow an interme- 
diary procedure. We compared the chi-square deviations of 
the data from the fitted curve, divided by the number of 
degrees of freedom (hereafter x 2 /dof). However, even when 
this parameter is smaller for a model with a break than for 
the straight line, there is still a probability of the straight 
line being the "real" one. We must then be able to attribute 
relative probabilities to the different functions that can be 
fitted. We performed Monte-Carlo simulations to answer this 
question. 

In the tests we made, a random noise was introduced 
on the top of a well behaved law, that we will call the "real" 
law, and then attempted to recover this original law by an 
automatic fitting procedure. For instance, it may happen 
that even starting from a real straight line, the fitting pro- 
cedure tells us that a third order polynomial is the best fit, 
based on the chi-square test. By repeating the experiment 
hundreds of times, it is possible to infer the probability of 
such an event to occur. 

The metallicity distributions were simulated starting 
with a straight line which represents the average oxygen 
abu ndance for spiral gala xies derived from the tables 1 and 
3 of IZaritskv etaO (|l994h (Eq-HJ. 

Z(r) = -0.09 -r + 8.98 (1) 

In a first part of the tests we generated data points 
with a random dispersion of 0.3 dex relative to this line. 
This dispersion is typical of the observed data. In a second 
part we created a new "real" radial metallicity distributions 
with a minimum with a depth of 0.3 dex with respect to 
(Eq. [TJ, by adding a cosine shape function in the middle 
of the distribution. The width of the minimum is similar 
to the typical uncertainty for the corotation radii (3 kpc). 
After that we generated data points with 0.3 dex dispersion 
with respect to the new distribution. The experiment was 
repeated using two connected lines changing the slope in 
the middle of the distribution, instead of a minimum, as 
being the original function. The change of slope was taken 
to occur at 0.3 dex below the straight line. 

For each original or real function we performed 200 sim- 
ulations by generating different sets of data points, and ana- 
lyzed these sets as follows. The simulated data were analyzed 
by automatically fitting different functions: a linear fit, a 4th 
order polynomial, or two connected straight lines with differ- 
ent slopes. For each fit we recorded the differences between 
the x 2 /dof of the fitted models and that of the single straight 
line fit. The probability of having a correct identification as a 




Figure 2. Results of Monte-Carlo simulations. A random distri- 
bution of points is generated around an original function (two 
straight lines for the upper panel, a function with a minimum 
for the lower panel and a straight line for the two panels). The 
original function and a straight line are fit to both set of data in 
each panel (with and without inflexions), giving us x 2 /dof(fit) and 
X 2 /dof(line) for each case. The histograms compare the frequency 
of the cases for which the correct function is recovered (repre- 
sented with dropping lines hatch pattern) and the frequency of 
the straight line being the best result (rising lines hatch pattern), 
both as a function of x 2 /dof(fit)-x 2 /dof(linc). Each histogram is 
the result of 200 simulations. The continuous line represents the 
probability of the correct function to be recovered, which is read 
in percentage in the scale. 

function of the difference of x 2 /dof are presented in Figure [2] 
Let us call a the parameter [x 2 /dof(model)- x 2 /dof( straight 
line)] 

As it could be expected there is a range of the param- 
eter a, corresponding to small differences of the chi-squares 
of the two fits, for which the quality of the fits is similar, 
and for which the probabilities of one or the other functions 
being the real one are similar. However, when a becomes 
larger than 0.025, the initial model function becomes the 
most probable one, compared to a straight line. For this rea- 
son, we will classify as "strong point" the radial metallicity 
distribution of galaxies for which a > 0.025. 

Interestingly, we found that even in situations where the 
above argument based on the value of a does not allow us 
to eliminate the possibility of a straight line, in most cases 
the polynomial fit provides a quite good result, in the sense 
that it recovers the correct position of the minimum, when 
the original function had one. 

4.2 Adopted procedure 

The radial abundance distribution data of each galaxy were 
treated by a fully numerical and automatic analysis proce- 
dure, as follows: 

(1) The observed metallicity distribution was fitted by sin- 
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gle straight line, two connected straight lines (4 degrees of 
freedom), two straight lines with a step (5 dregrees of free- 
dom), and then 4 polynomials, of order 2 to 5 (the number 
of degrees of freedom varies from 3 to 6). We verified that 
there is no need to go to higher orders. The uncertainties 
were taken to be equal for all the data points. 

(2) For each fit the xV dof 

was evaluated according to 
IPress et alj (|l992h - A descr iption of a similar a pplication 
of this recipe can be seen in Sc arano et al.l (|201lf ); 

(3) For the two connected straight lines and the two straight 
lines with a step, the point of connection or the step position 
respectively was assumed to be the point where the gradient 
break occurs, included in Table [2] For the polynomials we 
identified the minima and inflexions by studying the deriva- 
tives of the fitted curves. If the step function fitted the data 
better than the other fits, the inflexions points of different 
polynomials were used to determine the error on the radius 
of the metallicity break. If the two connected straight lines 
fitted the data better than the other fits, the polynomial 
minima were used to calculate the error on the radius of the 
metallicity break; 

(4) Taking into account only the results for which the 
X 2 /dof was lower than that obtained with the linear fit, we 
registered the radius of the minimum for each function fitted. 
If a fitted polynomial presented more than one minimum, 
we only considered the one situated near the minimum of 
the other polynomials. With these values we estimated the 
position of the break in the gradient and its dispersion. We 
adopted the weighted mean of the positions, using as weights 
the inverse of the xVdof. 

Figure [3] illustrates the application of this procedure to 
the galaxies NGC3359 and NGC5236. These galaxies are 
representative of the two classes of breaks found in the lit- 
erature: smooth breaks and step breaks. For NGC3359 we 
can see the general behaviour of the first class of objects: 
the radial metallicity distribution changes smoothly, and it 
is better fitted by two connected straight lines. Higher order 
polynomials also produce good fits, with minima close to 
the intersection of the two connected straight lines. Lower 
order polynomials also tend to reveal breaks, but depending 
on the data dispersion and the steepness of the slopes, the 
minima found are farther than those obtained with other 
models. On the other hand NGC5236 represents the second 
class of breaks. In general, for these cases, only the step 
function fits the data properly. As it could be expected, the 
minima of polynomial fits are not correlated with the step 
position. However, the inflexions of the polynomials tend 
to be near the position where the gradients of abundances 
changes abruptly. 



5 RESULTS 

We performed the analysis of the galaxies listed in Table [T] 
following the procedures described in the previous section. 
The data needed for this study were found in the literature, 
the references being given in Tables [ft]] and [ftj (Appendix) . 
The results are presented in Table [2] where the breaks or 
minima in the metallicity distribution obtained by the au- 
tomatic fitting method are given. 

The main result of the present work is the existence of a 
clear correlation between corotation radii and breaks in the 



metallicity distribution, shown in Figure [4] left panel. The 
linear fit to the data presented in this figure is very near the 
one-to-one match, which means that the two occurrences 
are spatially associated in the galaxies studied here. As ex- 
pected from the inverse relation between corotation radius 
and pattern speed, a good correlation can also be seen in 
a plot of angular speeds (angular velocity of the rotation 
curve at the radius of metallicity break versus angular pat- 
tern speed) , in the right panel of the figure. The uncertainty 
is large for some points in this case, among other reasons, 
due to the non-linearity in the process of conversion of the 
data from corotation radius to angular velocity, especially 
when no uncertainties were mentioned in the literature and 
we had to assume the median uncertainty in the corotation 
radii of the other references (sigma ~ 2kpc). 

In spite of our limited sample some side results can be 
seem at bottom of Figure [4] The distribution of pattern ro- 
tation speeds shows a broad peak at 24 km/s/kpc, which 
reveals that our Galaxy has a "typical" pattern angular ve- 
locity. This peak would be even more prominent if only the 
"strong" points were shown in the histogram. This result is 
certainly a clue to understand what determines the value of 
fi p in galaxies. On the other hand, the distribution of coro- 
tation radii (without any normalization to the size of the 
galaxy, in our work) is somewhat uniformly spread between 
about 3 kpc to 17 kpc. Considering that lar ge spiral galax- 
ies m ay reach typically a radius of 30 kpc (|Tullv fe Fisherl 
Il987l ) which should correspond to the outer limit for spiral 
arms (the Outer Lindblad Resonance), the cutoff seen in the 
corotation radius histogram of Figure [3] is coherent with the 
expectation that the corotation radius must be near midway 
between the Lindblad resonances. 



6 COMMENTS ON RESULTS FOR 
INDIVIDUAL GALAXIES 

We next present comments on the 27 galaxies included in 
Tables [T] and which are the galaxies for which we found 
determinations of both metallicity gradient and corotation 
radius. Some of these galaxies were not included in the cor- 
relation, for reasons that are explained in the notes, but they 
are included in our statistics on corotation radii. Of course, 
there are many more galaxies in the literature, not included 
in our study, for which only one of these two sets of data 
was found. One example is NGC1512, for which there is a 
nice set of measurements of O abundance by Bresolin et al. 
(2012), which shows the presence of a break in the slope 
of the distribution at 12-13 kpc. However, we did not find 
a determination of the corotation radius in the literature. 
We classified the galaxies that remained in the correlation 
shown in Figure [4] in two classes, the "strong" and "weak" 
points. These denominations do not have any direct relation 
with the error bars which affect the points. We consider that 
we have a strong point when the following conditions were 
fulfilled: 1-) more than 16 measurements of O abundance 
at differen t radii are available , satisfying the condition pre- 
scribed bv lDutil fc Rovl (120011 ): 2-) at least two different cal- 
ibration methods were used to determine the abundance and 
3-) the breaks (minima, breaks or steps) in the abundance 
distribution derived from the different calibration methods 
are in agreement. Galaxies for which there are between 10 
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Table 2. Corotation radii and breaks of metallicity gradients. The index n is used to represent the galaxy in the graphs. Hcr and Q p 
are respectively the median of the corotation radii and the spiral pattern speeds found in the literature, discarding bad values, Rdz is 
the metallicity break obtained by the procedure written in the text, and Q^z is the angular velocity of the rotation curve at R^z- The 
"Notes" column present comments on the quality of the data or of the fit to obtain the metallicity breaks (S= strong point, W= Weak 
point, FP= too few points in the metallicity distribution, see text for explanation of S and W). The last column gives the references for 
the radial metallicity distributions according to Table fATl in the Appendix. 
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Galaxy 


Res 




RdZ 


^dZ 


Notes 


References 






[kpc] 


[km/s/kpc] 


[kpc] 


[km/s/kpc] 






1 


IC0342 


10.2 ± 3.1 


19.0 ± 7.4 


- 


- 


FP 


[76], [89], [125] 


2 


NGC0224 


15.6 ± 1.1 


18.1 ± 1.1 


15.3 ± 3.6 


18.2 ± 4.6 


W 


[10], [62], [89], [125], [133], [135] 


3 


NGC0598 


4.3 ± 1.4 


24.0 ± 6.1 


5.4 ± 1.3 


19.1 ± 4.4 


W 


[10], [22], [54], [62], [71], [76], [89], 
















[1091 f125l [1351 


4 


NGC0628 


4.6 ± 1.2 


19.2 ± 7.6 


4.8 ± 1.8 


18.4 ± 6.5 


w 


si loi ri7i f38i \76] rsoi rssi 
















[89], [109], [121], [122], [125], [135], [137] 


5 


NGC1068 


3.2 ± 0.1 


67.5 ± 2.1 






FP 


[36], [83], [89], [121], [122], [135] 


6 


NGC1232 


12.6 ± 3.0 


25.6 ± 6.1 


15.0 ± 2.1 


21.6 ± 3.1 


S 


[9], [17], [88], [89], [121], [122] 


7 


NGC1365 


20.7 ± 4.0 


11.5 ± 5.6 


20.4 ± 5.3 


11.6 ± 3.8 


s 


[9], [28], [62], [87], [89], [98], [125], [135] 


8 


NGC1530 


13.4 ± 5.0 


8.2 ± 5.0 


- 


- 


FP 


[72] 


g 


NGC1566 


8.8 ± 1.3 


23.6 ± 5.0 






FP 


[62], [109], [125], [135] 


10 


NGC2403 


5.1 ± 3.0 


25 ± 15 


6.6 ± 0.6 


19.2 ± 1.7 


W 


n 01 [401 [471 [761 [801 f87l 
L-luj, [iuj, L' U J) l ou J' L° ' J ; 
















[89], [109], [121], [122], [125], [135] 


11 


NGC2543 


6.0 ± 2.2 


25.9 ± 9.3 






FP 


[72] 


12 


NGC2903 


4.1 ± 0.1 


69.9 ± 2.7 


5.0 ± 0.8 


56.4 ± 10.5 


S 


[62], [76], [88], [89], [109], [121], [122], 
















[125], [135] 


13 


NGC3031 


11.0 ± 2.0 


18.2 ± 5.0 


9.8 ± 3.4 


21.7 ± 12.1 


s 


[10], [54], [62], [85], [89], [114], [125], [135] 


11 


NGC3319 


7.4 ± 2.0 


9.8 ± 0.7 


6.6 ± 0.7 


10.1 ± 0.3 


w 


[62], [135] 


15 


NGC3359 


6.1 ± 0.6 


20.7 ± 2.7 


7.9 ± 0.5 


16.2 ± 0.8 


s 


f73l [1341 


16 


NGC3621 


6.2 ± 2.0 


22.0 ± 4.7 


8.7 ± 1.0 


16.1 ± 1.8 


s 


[12] 


17 


NGC4254 


13.6 ± 2.2 


10.5 ± 1.6 


11.5 ± 3.0 


12.6 ± 2.2 


w 


f62l [761 [801 [881 [891 [1091 [1121 
L J' L J' L J' L J' L J' L J' L J' 
















[1 251 [1 351 


18 


NGC4321 


10.5 ± 1.3 


25.8 ± 3.0 






FP 


[28], [62], [76], [80], [88], [89], [109], 
















[112], [125], [135] 


19 


NGC4736 


3.1 ± 0.4 


56.0 ± 5.0 






FP 


[10], [62], [75], [76], [83], [89], [125], [135] 


20 


NGC4826 


1.2 ± 0.1 


56.2 ± 5.9 






FP 


[83] 


21 


NGC5033 


10.4 ± 0.7 


25.1 ± 1.0 






FP 


[62], [80], [89], [135] 


22 


NGC5055 


14.7 ± 4.4 


14.6 ± 3.8 






FP 


[62], [76], [89], [125], [135] 


23 


NGC5194 


5.6 ± 0.6 


45.6 ± 4.3 


6.8 ± 0.5 


37.0 ± 2.7 


S 


[8], [10], [62], [76], [80], [88], [89], 
















[125], [135] 


21 


NGC5236 


8.1 ± 2.0 


22.5 ± 6.5 


8.4 ± 3.0 


21.8 ± 9.9 


S 


[9], [10], [11], [62], [88], [89], [125], [135] 


25 


NGC5457 


16.7 ± 2.4 


5.5 ± 1.2 


17.7 ± 2.2 


5.9 ± 0.5 


s 


[28], [54], [61], [76], [87], [88], [89], 
















[109], [121], [122], [125], [135] 


26 


NGC6946 


4.2 ± 0.4 


53.8 ± 1.9 






FP 


[5], [38], [62], [76], [88], [89], [109], 
















[125], [135] 


27 


NGC7479 


14.6 ± 2.7 


12.8 ± 2.4 


16.6 ± 0.6 


11.0 ± 0.5 


S 


[72], [74], [89] 
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Figure 3. Radial oxygen distribution (top panel) and metallicity gradients (bottom panel) for the galaxies NGC3359 (left) and NGC5236 
(right). The models fitted to the data are represented by different line types, as indicated in the legend. 



and 16 measurements at different radii or only studied using 
one single calibration method, even showing a clear break in 
the abundance distribution, are considered as weak points. 
A few galaxies which fulfil the "strong" conditions were con- 
verted to weak points for some other particularity like non- 
even distribution of the points, inhomogeneity of the metal- 
licity sample, or other, as explained. 

IC0342: This galaxy was not included because only five H n 
regions were observed. 

NGC0224 (MSI): It is difficult to see the spiral structure 
of this galaxy, due to its high inclination. There are many 
abund ance measurements of HII regions. One set of data is 
that of lPilvuein et alJ l|2004h . However, with these data, the 
best fit obtained is linear (without breaks in the slope). We 
attribute this to the scarcity of data beyond 10 kpc, where 
the break occurs. A change of slope around 15 kpc can be 
seen i n the data of iBlair et alJ (1 19821 ) and of IWorthev et all 
l|2005l ). The data of Worthey (34 objects) include, in addi- 
tion to HII regions, stars observed with the Hubble Space 
Telescope, as distant as 50 kpc. The data are well fitted 
with two lines with a break at 15.2 kpc. The extended re- 
gion with very little change in metallicity, from about 15 to 
50 kpc, is well fitted by a straight line and makes the poly- 
nomial fits to look unsatisfactory. Although there is good 
agreement between the break and the corotation radius, we 
will consider this galaxy as a "weak" point, due to the in- 
clusion of different types of objects in the determination of 
the gradient. 

NGC0598 (M33): There are many discrepan t determina- 
tions of corotation for this galaxy. iNewtonl |l980) con- 



siders that the corotati on radius can be 1.2 or 7.0 kpc 
(scaled to our system). iPueraril (|l993h determines a pat- 
tern speed which corresponds to corotation radius of 
4.3 kpc. A rece nt work on the distribution o f plan- 
etary nebula by Magrini, Stanghcll ini fe Villaverl (2009) 
shows a gap at 3.8 kpc (also scaled), with a constant 
or slig htly rising oxygen abundance beyond that radius. 
Larger ([Vila-Costas fc Edmunds! (Il992h 5.9 kp c) and smaller 
|Elmegreen. Elmegreen fc Montenegro! (|l992l ) 2.9 kpc) radii 
can be found as well in the literature. We adopted the me- 
dian corotation radius of the data from the literature. This 
result is coherent with extension of the spiral structure. Con- 
cerning the metallicity distr i butio n, recent results h ave been 
produ ced bv lMagrini et all (|2007l ) (28 points) and iBresolinl 
|201lf) (46 points). In the work of Bresolin et al., there are 
some "best quality" points, but they do not reach the ra- 
dius of corotation; the other points, collected from the liter- 
ature, show a larger scattering and are best fitted by a single 
straight line. The data of Magrini et al. on planetary nebu- 
lae is best fitted by two lines with a break in the slope. This 
result is consistent with the break at the same radius, but 
considering that this result was obtained only for one cali- 
bration method, this galaxy is considered as a weak point. 

NGC0628 (M74): The metallicity distribution of this 
galaxy was studied using two sets of data inde- 
pendently. One set is that of using the data by 
iRosales-Ortega et all (|201lh . who obtained abun- 
dances of 96 regions using integral field spectroscopy 
with 4 different calibration methods. If we consider 
only the new data obtained by these authors, for 
two of these calibrations the fit of two lines with 
a step reveal a very smooth break in the metal- 
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Figure 4. Correlations between the corotation radii and the metallicity breaks radii (left panel) and the spiral pattern speed and the 
angular velocity expected at the metallicity breaks radii (right panel). We labelled the data points according the index sequence of 
galaxies in Table [2] ("n" column). Our Galaxy is represented by a "G" letter; the dashed line corresponds to the perfect match line, 
the dashed-dotted line is the linear fit considering all dataset; the continuous lines take into account only galaxies with the strongest 
evidences of breaks in their oxygen distribution (see text). The histograms at the bottom part of each correlation graph show the number 
of objects per bin of data (corotation or spiral pattern speed), taking into account all galaxies presented in Table 3 (bin size equals the 
mean of uncertainties of these data). 



licity gradient near 3.6 kpc. Better results are ob- 
tained when data from the literature are added 
to the new observations, as proposed by Rosales- 
Ortega et al. The best fitted model to all calibration 
methods present a step near 6.1 kpc. However, the 
histogram of radial breaks found as a function of 
the galactrocentric distance exhibit two peaks: one 
near 5 kpc and another near 8.5 kpc for all cali- 
b rations. The oth e r set of data that we analysed is that 
of iMoustakas et al.l (|2010l ). which present two distinct cal- 
ibration methods, with 34 p oints. These are the t h eoret- 
ical calibration published by iKobulnickv fe Kewlevl (|2004l . 
hereafter KK04) and th e empirical calibration method by 
IPilvugin fc Thuanl (|2004 hereafter PT05). For the first cal- 
ibration the data is best fitted with two straight lines or two 
lines with a step, with the breaks at 6.3 and 6.9 kpc respec- 
tively. The results of the PT05 calibration provide better 
fits than those of the other calibration, and are best ad- 
justed with polynomials of order 4 and 5, followed by 2 
lines, and two lines with a step. The average for this 
calibration, weighted by \ 2 > i s a break at 4.8 kpc. 

This radius agrees with the corotation determina- 

I 1 } TJ\ I 1 

tion of Cepa fc Beckman (1990), Elmegreen et al. 

(1992) and IVila-Costas fc Edmunds! (|l992l ) which 



give an average of 4.6 kpc. Interestingly, for the 
data of this calibration the fit of two lines with a 
step indicates the existence of a step at about 7 
kpc. This result is not much different from the sec- 
ond peak in the histo gram that we mentioned fo r 
the data analyzed by Rosale s-Ortega et al.l (|2011). 
Taki ng into account t he ro tatio n curves provided by 
iKamphuis fc Briggsl (| 19921 ) and iDaigle et all (|2006T ) 
we see that the corotation at 4.6 kpc provides a spi- 
ral pattern speed of 19.2 km/s/kpc. In this case the 
outer Lindblad resonance (the limit where the spiral 
arm can exist) would be placed at 7.2 kpc, very near 
the end of the visible spiral arms. Since beyond the 
end of the spiral structure there would not be an ef- 
ficient and continuous process to enrich the medium, 
a break in metallicity is also expected at the outer 
Lindblad resonance, and we could attribute to this 
resonance the outer break seen in the histogram of 
radial breaks. However we do not exclude the pos- 
sibility of corotation being between 7.5 and 8.5 kpc, 
since this value would be more consistent with the 
diffuse e xtension of the sp iral structure seen in the 
H I line (Wa lter et al.ll2008h . Also we cannot discard 
the superposition of multiple pattern speeds. Due to 
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the lack of close agreement between the calibration methods, 
we classify this galaxy as a weak point. 

NGC1068: Only 9 H II regions have been measured, and 
in addition, they are poorly distributed along the disk. We 
exclude this galaxy from our correlation analysis. 

NGC1232: G al axy sampled by 16 H 1 1 regions by 
Pilyugin et al l (|2004l ), originally fr om Ivan Zee (1998), to 
which new sources were added by iBresolin et all (|2005h . 
Since there are points in common in the two samples, one 
can see that the two calibrations are in good agreement, 
so that we can add the samples, resulting in a total of 23 
H II regions. The data is well fitted by polynomials of order 
3,4 and 5, and by two lines. The radius of the break from 
these different fits were averaged. This galaxy can be con- 
sidered as a strong point; the coincidence of the break with 
the corotation radius is satisfactory. 

NGC1365: Pilyugin collected about 70 points in the lit- 
erature for NGC13 65, most of which (55) are due to 
iRov fc Walshl (|1997| ). who were the first to detect the break 
in the metallicity slope of this galaxy. Later [Pilvuginl (|2003l ) 
put in doubt the reality of the break, which could be an 
artefact due to the existence of two branches in the calibra- 
tion of the R23 method. We have several arguments telling 
us that the break is real: 1-) Roy and Walsh compared 3 
different methods of calibration for his data, and not all of 
them necessar il y affe cted by the same systematic error; 2-) 
IBresolin et al.l (2005) added new data points which allow a 
comparison of the P-method with the Te-method. The two 
methods give the same values of 12+ log[0/H], between 8.6 
to 8.4 dex, in the range of distances from 15 to 20 kpc. This 
suggests that the abundances found beyond 20 kpc, which 
are in the same abundance range, were not affected by a 
systematic error in the Pilyugin's data; 3-) In the same pa- 
per in which the reality of the break is discussed, it is said 
that problem of calibration could affect abundances smaller 
than 8.2 dex. However, almost all the points are above this 
limit, in the case of NGC1365, even taking in to account 
the data dispersion of 0.2 dex, so there is no reason to be- 
lieve that all data points could be affected by calibration 
errors. Consequently, we will use not only the data of Pi- 
lyugin, presented in his Figur e 2b, for our analysi s, but also 
all calibrations presented bv lRov fc Walshl ij 19971 ). One dif- 
ficulty that one has to face is that there are two solutions. 
The data is well fitted by polynomials of 3rd, 4th and 5th 
order, and by 2 lines, which give a minimum at about 20 
kpc. This is the minimum which is shown with 2 lines in 
Pilyugin's Figure 2b. However, the data is also marginally 
well fitted with two lines with a step at 16.6 kpc, close to 
inflexion points revealed by the same 3rd, 4th and 5th order 
polynomials. It should be remembered that in the presence 
of a step, we must look for the inflexion points of the poly- 
nomials, and not for the minima, as explained in Section 
2. There seems to be a robust case for the presence of the 
step, more clearly seen in the data based on R23. Concerning 
corotation, there are several measurements in the literature, 
the median being 13.9 kpc, but most of these values are in- 
compatible with the extension of the sp iral structure if we 
take into account the rotation curve by ISofud (|l996l ). The 
only value compatible with the observed spiral arms is an 



alternetive value mentioned bv lVera-Villamizar et al.l (|200lh 
of 2 0.7 kpc (re-scaled ) . A re cent detailed study of NGC1365 
By lElmegreen et all (|2009l ) is of particular interest, since 
they detected a radial flow of gas along the spiral arms, in 
opposite directions on the two sides of corotation. The large 
error bar attributed to the radius of the break accounts for 
the two solutions (about 16 kpc and about 20 kpc). 

NGC1530: Its metallicity distribution has only three ob- 
served H II regions, so that it was discarded from our anal- 
ysis. 

NGC1566: Grand-design galaxy with six H II regions ob- 
served, so that it was discarded. 

N GC2403: Flocculen t galax y with 45 H II regions compiled 
by iMoustakas et all (|2010l ). For one of their calibrations 
(KK04 method) all models proposed here fit the data better 
than a straight line and show the presence of a break, while 
for their other calibration (PT05 method) the data is best 
fitted by a straight line. For this reason we classify it as a 
weak point. There is only one reference giving the cor otation 
radius of this galaxy (|Vila-Costas fc Edmundsll 19921 ). which 
is compatible with the radius of the break in metallicity that 
we found. 

NGC2543: Galaxy with too few observed H II regions. Al- 
though it was discarded from our analysis, it is interest- 
ing to mention that at least visually there is a step in the 
metall icity distribution of this galaxy at the sam e position 
where rVera-Villamizar. Puerari fc Dottorll ((2003) estimated 
the corotation radius position using Fourier analysis. 

NGC2903: There are two evaluations of the corotation ra- 
dius for this galaxy, coincident with each other and with the 
metallicity break found at about 5 kpc using the polynomi- 
als of order 3, 4 and 5, in addition to the two connected lines. 
The metallicity distributio n is well sampled by 36 H II re- 
gions l|Pilvugin et al.ll2002l ) and the results seem to be com- 
patible with the spir al structure extension . The same result 
is obtained using the lZaritskv et all (|l994h . which make this 
galaxy a strong point in the correlation presented in Figure 
[4] We cannot exclude the possibility of the corotation be- 
ing at 9 kpc, where our polynomial fits presents an inflex- 
ion, but there are no references in the literature supporting 
that. This would put the corotation radius at the extremity 
of two well behaved symmetrical spiral arms, like it happens 
in other galaxies, and would bring the pattern speed to 31 
km/s/kpc, a more usual one than 57 km/s/kpe. 

NGC3031 (M81 ): The estimates of the corotation radius for 
this galaxy is sometimes mentioned to be around 8 kpc and 
other times around 1 4 kpc, as can be seen in the work by 
iKendall et al.l |2008al ). recalibrated according the projection 
parameters assumed in this work. In spite of this, based on 
its frequency diagram, the extension of its spiral structure 
favors a spiral pattern speed of 20 km/s/kpc. This speed re- 
sults in a corotation radius near 8 kpc, compatible with the 
breaks found using the polynomial fits of order 3 to 5 and 
with the tw o straight line fits of the radial abundances de- 
termined bv lStanghellini et al.l (|2010l ) using planetary neb- 
ulae. These results are also consistent with the breaks found 
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using the radial abundance distribution bv lPatterson et all 
(2012) for H n regions with the calibrator KK04, given 
in their Figure 9, which also includes data from other 
authors. For this reason this galaxy is assumed as a strong 
point in our correlation. 

NGC3319: This is another example of a highly inclined 
galaxy. In contrast with NGC0224 it has well defined break 
identified by all models that we fitted to its metallicity dis- 
tribution. Unfortunately it is sampled only by 13 H II re- 
gions for which only one calibration method was applied 
l|Zaritskv et al|[l994l ). so it is represented by a weak point 
in the correlation in Figure [4] Its rotation curve is almost 
a rigid body rotation curve (accordin g the results in optic 
bvlGarrido et al.ll2002l and in radio bv lMoore fc Gottesmanl 
Il99sl ). Since its spiral pattern speed and the angular ve- 
locity associated to the rotation curve are similar along a 
significant part of the disk it may explain the "ring" with 
less H II regions in this galaxy. Note that the spiral pat- 
tern speed mentioned in the literature is compatible with 
the spiral structure extension. 

NGC3359: Galaxy with 94 H u region s obs erved by 
iMartin fc Rovl (| 19951 ) and lZahid fc Bresolinl <|201ll ). In both 
references the authors found breaks in the radial metal- 
licity distribution using different calibration methods and 
associated it to a possible effect of the corotation radius. 
IZahid fc Bresolinl confirmed the statistical significance of the 
flattening using the an F-test. Similarly, according our pro- 
cedure, polynomials of order 4 and 5, two connected lines 
and two straight lines with a step fit the data better than 
the linear model. For all this cases a break in the metallic- 
ity distribution is fou nd around 7.9 kpc, near the corotation 
radius determined by iRozasl (120081 ). For these reasons this 
galaxy is represented as a strong point in the correlation 
shown in Figure U 

NGC3621: The radial o xygen abundance of this galaxy was 
extens ively studied by iBresolin. Kennicutt fc Ryan- Weber! 
i|2012l) using different calibrators for its 73 H II regions. 
These authors found a significant flattening in the metallic- 
ity gra dient of this galaxy after R25 using the same proce- 
dure of lZahid fc Bresolinl i|201ll ). Our procedure a grees with 
their results and two connected lines and two straight lines 
with a step fit the data and exhibit an evident break near 8.7 
kpc. These results from different analyses turn this galaxy a 
strong point in the correlation shown in Figure 4. Unfortu- 
nately there a re not many studies on the corotation radius 
of this galaxy. iTamburro et al.l l|2008l ) determined its spiral 
pattern speed by deriving geometrically the timescale of the 
star formation driven by its spiral arms. For this purpose 
they used SINGS and THINGS data sets to verify the angu- 
lar offsets between the 24 /im and the H I emission tracers. 
The original result of this work points to a corotation radius 
of 2.2 kpc. However these authors tag this value as a result 
of a bad quality fit. We derive a better value (6.2 kpc) by 
discarding the most d iscrepant data from the analysis by 
ITamburro et al.l (|2008l ) an d comparing the resu lts with the 
observed rotation curve bv lde Blok et all (|2008l ). 

NGC4254 (M99): Galaxy sampled by 18 H 11 regions with 
oxygen abundance determination listed bv lMoustakas et all 



(2010). The objects are relatively well distributed along the 
disk but the observations do not go too far from the me- 
dian corotation radius mentioned in the literature. None 
of these references on the corotation fully agree with each 
other, but the metallicity break found using a 3rd order 
polynomial and two lines with a step on the data cali- 
brated with th e PT05 method a r e compatible with the me - 
dian results bvlElmegreen et all (|l992l ). iKranz et all (|200ll ) 
and iGonzalez fc Grahaml (Il996l ). which we adopted in our 
correlation. On the other hand, the break radius found 
using the KK04 calibration and two connected lines fit 
is consistent with the cor otation radius determined by 
IGonzalez fc Grahaml (|l996l ) and we cannot discard their 
value as the real position of the corotation radius. Since the 
results from different calibrations do agree with the presence 
of a break, but they do not agree with a specific position in 
the radial oxygen distribution, this galaxy is represented as 
a weak point in our correlation. 

NGC4321 (M100): This is the most explored galaxy regard- 
ing the evaluation of the corotation radius and the spiral pat- 
tern speed. A histogram of the corotation radii taken from 
the literature reveals a double peaked distribution, indicat- 
ing that two spiral pattern modes could be superimposed. 
The extension of the spiral structure is compatible with the 
corotation radius of the external peak, which is also con- 
sistent with the metallicity breaks found with 4th and 5th 
order polynomials. Unfortunately the abundance distribu- 
tion of this galaxy is sampled by only 10 H 11 regions, so we 
had to discard it from our analysis. 

NGC4736 (M94): Galaxy with peculiar arms very extended 
in relation the brighter part of the disk where a ring can be 
observed. The sample of H 11 regions consist in only 8 ob- 
jects, all of them very concentrated in the internal regions of 
the galaxy, and for this reason this galaxy was not included 
in our analysis. 

NGC4826: Another galaxy with few H 11 regions (8) studied 
along the galactic disk. The objects are very concentrated 
at galactic radii internal to the expected position for the 
corotation radius. Therefore this object was also excluded 
of our analysis. 

NGC5033: This galaxy presents an asymmetric spiral struc- 
ture with distinct pitch angles. Considering that there is 
only 8 H 11 regions observed along this galaxy we did not 
take it into account in our correlation analysis. 

NGC5055 (M63): Flocculent spiral galaxy with only five H 11 
regions studied in terms of metallicity. Hence, it was not pos- 
sible to perform a reliable study of its metallicity distribu- 
tion, and consequentely we discarded it from our correlation 
study. 

NGC5194 (M51): Grand-design spiral galaxy with 16 stud- 
ies relative to its corotation radius. The values are bimodally 
distributed around 5.5 and 11 kpc (near the outer limit of the 
spiral structure) . There i s a very precise study on the metal- 
licity distribution by Bre solin. Garnett fc Kennicutt Jrl 
(2004), using temperature-sensitive line m ethods in 10 H 11 
regions and also an extensive study by iMoustakas et al] 
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l|201Cft . using the KK04 and PT05 methods in 32 H II regions 
compiled in the literature. In spite of the small amount of 
H n regions in the first case, breaks in the metallicity dis- 
tribution were found using polynomials of order 3 and 4, 
with values that are coincident with the median corotation 
radii from the literature. Breaks at the same position were 
found in the second case using fits of 3rd order polynomial, 
two connected lines and two straight lines w i th a step. Using 
the sa me method of iMishurov et ail (|2002l ). lAcharova et al.l 
(2008) determined the spiral pattern speed, of which the cor- 
responding corotation radius is coincident with the radius 
of the break. For these reasons this galaxy is considered a 
strong point in the correlation shown in Figure [4] 

NGC5236 (M83): This g alaxy is possib l y one of the best 
cases of the present study. iBresolin et al.l (|2009f ) obtained O 
abundances for 49 H II regions, with low scattering, which 
combined with data from 23 H II regions previously ob- 
served shows a clear change of slope associated to a step 
at 5.5 arcmin (8 kpc). There are very discrepant determi- 
nations of the corotation radius for this object. To make a 
choice, we must first examine some aspects of the morphol- 
ogy of this galaxy. 1-) The rotation curve is fiat to large 
distances, which indicates that the galaxy is not disturbed 
by any strong external i nfluen ce; 2-) A seminal paper of 
lElmegreen fc Elmegreenl (|l995l , hereafter E&E) shows that 
almost all spiral galaxies have a well defined 2-arms struc- 
ture up to a distance of 0.5 R25. We see the end of this struc- 
ture at about 3 arcmin from the centre (also determined by 
E&E for this galaxy), which places R25 and the corotation 
radius at about 6 a rcmin from the center. 3-) A map of 
the gas distribution l|Crosthwaite et al"1l2002r i shows a ring 
void of gas with radius 5.5 arcmin, which resembles the one 
which exists at coro tation, in our Galaxy. As discussed by 
lAmores et all 12009). such a gap is not only observed in our 
Galaxy but also predicted by theory as an effect of corota- 
tion. In general the corotation radius mentioned in the lit- 
erature is relatively smaller than th e indicators mentio ned 
here. A same grou p of researchers l|Talbot et al . 1979 and 
Ijensen et al.lll98lh obtained two discrepant values of fi p (22 
and 53 km/s/kpc, recalibrated accoriding our projection pa- 
rameters) using different methods. In their 1981 paper, these 
authors analyse the reasons for this discrepancy and state 
their preference for the largest value (which implies a smaller 
corotation radius). They comment that the first method re- 
lies on the fact that the star formation rate is proportional to 
(Q — fl p ) and reject this idea, although it gives a convincing 
fit of the data. Note that this proportionality is precisely the 
one we verify in our Galaxy. On the other hand, the argu- 
ment that they present in favour of the largest Q p is the ab- 
sence of a significant drift of the young stars with respect to 
the arms, near the inner end of the arms. They consider that 
the drift should be proportional to (Q — £l p ); if the drift is 
very small, the observed region must be close to corotation. 
However t hey do not take i nto a ccount the mechanism de- 
scribed bv lLord fc Kennevl l|l99ll) for M83 in which the gas 
gets trapped in the bottom of the channel-shaped potential 
traced by the arms, and flows along them. Consequently, the 
newborn stars have initial velocities along the arms, so no 
drift would be observed. Therefore we adopt 22.5 km/s/kpc 
as the correct solution. This places the corotation at 8.1 kpc, 
in agreement with our previous arguments. 



NGC5457 (M1 01): This g alaxy has a well sampled metallic- 
ity distribution IPilvuginl (|2003r i with 76 H II regions, which 
presents a clear flattening in the external regions. All the 
models that we use (polynomials, 2 lines, or two lines with 
a step) give a better fit than a single straight line, and point 
to a break at 19 kpc. However, in his 2003 paper, Pilyu- 
gin questions the existence of the bends in the slopes of 
the oxygen abundance in spiral galaxies. The best data for 
NGC5457, analysed with Pilyugin's calibration method, is 
shown in his figure lb, where we can see that all the points 
have abundances larger than 12 + log(0/H) = 8.2. In the 
same paper it is said that Pilyugin's formula is valid above 
this limit (8.2). If we take into account the uncertainties 
and the data dispersion (0.2 dex) we see that in spite of 
many data points be ing close to the confid ence limit they 
cannot be discarded. iKennicutt et alj (|2003T ) presents a set 
of measurement of 20 points, determined using the "direct" 
calibration, for which the radial oxygen distribution is sig- 
nificantly less flat than in the previous case. For these data, 
fitting with 2 lines, or with 2 lines with a step, the results are 
better than with a single line, and give breaks at 17.7 and 
17.9 kpc respectively. The median value of these breaks are 
in relative agr eement with th e corot ation radius found using 
the results bv lRoberts et all l|l975h . which is also compati- 
ble with the extension of the spiral structure. This last argu- 
me nt forces us discard th e val ues of the spi r al pat tern speed 
by iMeidt fc Rand! (j2005T l and ISheth etTal] i|2002ft . Since we 
found two sets of data in relative agreement and based on 
different calibration methods, this galaxy is considered as a 
strong point. 

N GC6946: Galaxy s a mpled by only 7 H II regions compiled 
by iMoustakas et all (|2010l ) and for this reason it will not 
included in our correlation analysis. 

NGC7479: This object is also one of the most extensively 
studied in term of its corotation radius, with almost 20 ref- 
erences on it. Most of the corotation radii are lower than 11 
kpc, all of them incompatible w ith the spiral structure ex- 
tensio n, exc ept for the r e sults b v lCourtes fc Dubout-Crillonl 
|l97ll) and iFathi et al] (|2009l ). In terms of radial oxygen 
abundance the data is almost excl u sively from the obser- 
vations by iMartin. Lelievre fc Rovl (|2000l ) of 69 H II re- 
gions candidates, from which at least 63 were suitable for 
abundances studies with four different calibrations meth- 
ods. They all look similar, which suggests that the step seen 
in the radial distribution is not an effect of a calibration er- 
ror. Analyzing the distribution with the smallest scattering 
in the data, we find that two lines with a step at 17 kpc 
and two connected lines with break at 16.2 kpc give better 
fits than a single line. Considering that the mean metallicity 
break radii obtained from different calibrations methods are 
similar, this galaxy is considered as a strong point in the 
correlation of Figure [4] 



7 DISCUSSION 

The existence of a minimum of metallicity at corotation, and 
the corresponding rising metallicity in the external parts of 
a number of spiral galaxies, can be understood under one of 
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the simplest assumptions that can be made, namely that the 
star formation rate (SFR) i s proportional t o [iR \ Q — Sl p |, 
or similar expressions, eg. IMishurov etaH l|2002l ). This is 
the surface density of gas /i times the velocity at which the 
gas penetrates the spiral arms, or in other words, the rate 
at which the gas feeds the star-formation mac hine. A simi- 
lar rec ipe for star formation was proposed bv lWvse fe Siikl 
(1989), but they did not give importance to the fact that 
one must use the modulus of the frequency difference. The 
relative velocity of the gas of the disk with respect to the 
spiral pattern is zero at corotation, and reverses its direc- 
tion at larger radii. The oxygen abundance at a given radius 
is proportional to the integral of star-formation rate with 
time, since O is mostly contributed by type II supernovae. 
These stars are short lived and do not have time to migrate; 
almost as soon as they are born they deliver new oxygen into 
the interstellar medium. Beyond the corotation the side of 
the arms on which the gas penetrates is reversed, but the 
modulus | Q — Q p | increases again, and consequently also 
the rate of enrichment of the medium in oxygen. 

This idea of proportionality between the amount of gas 
which enters the arms and the number of stars formed is 
specially easy to acce pt if we a dopt the classical descrip- 
tion of spiral arms of iKalnaisI (|l973l ). in which the arms 
are due to crowding of stellar orbits in some regions of the 
galactic disk. Therefore, the arms are "massive" entities, and 
behave like elongated potential wells or grooves in the ax- 
isymmetric potential. The interstellar gas gets trapped and 
accumulates into the elongated wells until the Jeans mass is 
reached and star formation starts. In this process the flux 
of gas reaching an arm is in some way "integrated" over the 
time. The above explanation gives us a rough picture of pro- 
cesses which might be going on. The description is a little 
oversimplified, however, since it is not only the velocity of 
galactic rotation relative to the spiral pattern that forces the 
gas to penetrate into the arms. The gravity of the arms, by 
itself, is able to attract the gas. 

The presence of a step in metallicity at corotation, in 
addition to the change of slope, as we observed in several 
galaxies, requires another explanation. In our view, this is 
a consequence of the fact that the two sides of corotation 
are isolated one from the other, with the radial gas flow 
induced by the spiral pattern going inwards inside corota- 
tion and outwards on the outer side. As a consequence, the 
two sides have no exchange of gas and evolve independently. 
This interpretation of the step was discussed by L+8 in the 
case our Galaxy. Such opposit e flows are not only predicted 
by th e theory of spiral arms IMishurov. Amores fc Lepind 
( 2009 ), but confir med by hydrodynamic simulations (eg 
ILepine et alj (|200ll ). Figure |4l) and directly observed , like i n 
the case of NGC1365 [Elmegreen. Galliano fc Alloinl |2009h . 
In addition to be isolated one from the other, the explana- 
tion requires that the two regions have different SFR. This 
is easy to accept since we always see the inner regions of 
galaxies having higher stellar and gas densities. 

Both the positive gradients in external regions of spiral 
galaxies, and the metallicity steps at corotation radius, of 
which we confirm the existence in this work, have an impor- 
tant consequence for the understanding of the evolution of 
spiral structures. A relatively long time is required to build 
up these features. This is a strong argument against the con- 
cept of "transient" spiral arms. If the corotation radius did 



not remain for a long time at a same radius, the disconti- 
nuities in chemical abundance would be smoothed out. L+8 
estimated the minimum lifetime of the present spiral struc- 
ture of our Galaxy to be about 3 billion years. This is also a 
typical number for the galaxies investigated here, based on 
the same arguments. 

The correlations presented in Figures 2] in some way 
validate the methods used to find the corotation radius 
in the literature, and the method that we used to locate 
breaks in the metallicity gradients. If any one of these meth- 
ods failed frequently, the correlations would be much worse. 
This is also an argument against the idea that galaxies of- 
ten have multiple pattern speeds. If this were so, we could 
have the corotation radius of one of the pattern speeds con- 
fronted with the metallicity features of another one, and 
this would destroy the correlation. It seems that at least a 
long lived dominant pattern speed exists for each galaxy. 
On the other hand we cannot discard the possibility of some 
mis-association between breaks in the radial metallicity dis- 
tribution and the other resonances, specially the Lindblad 
resonances and the 4:1 resonance. 

It is worthwhile, at this point, to pay some attention to 
minor deviations with respect to the one-to one relation- 
ship. One may note a tendency of the metallicity break 
being at a slightly (10 %) larger radius than corotation, 
in Figure |4^i. In agreement with this tendency, the fitted 
straight lines to the angular speeds (Figure [4p) have a 
slight bias in the sense of spiral pattern speeds being faster 
than the one-to-one match line. The minor deviations of 
the two correlations could be attributed, for instance, to 
a systematic offset of the metallicity break towards a ra- 
dius larger than corotation. The position of the minimum 
of the star formation rate (which originates the minimum 
of metallicity that we detect) can be influenced by the gra- 
dient of gas density. One can note for instance in the re- 
sult of a hydrodynamic simulation of the effect of the spiral 
potential perturbat ion on the gas density, in Figure 4 of 
ILepine et all |200ll) paper, that the resulting minimum of 
gas density is at a slightly larger radius than corotation. On 
the other hand, with respect the spiral pattern speed, the 
deviations with respect to perfect match could be attributed 
to systematic errors in the determination of the corota- 
tion radii, specially when th e Tremaine- Weinberg method 
IITremaine fc Weinberg! [l984l TW-method hereafter) is ap- 
plied. Apparently the TW-method provides pattern speeds 
systematically higher than the other methods, especially 
when it is applied to bars. If the spiral arms and the bar 
have the same pattern speed, as suggested by the clear 
bar-arm connection in mos t gal axies (see discuss ions in 
iKormendv fc Normanl Il979l and iBlock et all |2004| ). there 
could be a source of errors associated to the TW-method 
which results in higher velocities as seem for the objects 
with higher spiral pattern speed in [4] Possibly the TW- 
method is affected by the dispersion of the stellar veloci- 
ties. This suspicion is supported by the results of an anal- 
ysis of the corotation radius scaled by the semi- major axis 
obtained bv lCorsini. Aguerri fc Debattistal (|2004l ) using the 
TW-method, compared to the central velocity dispersion 
compiled from the hyperledaQ database. Figure [5] shows 
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Figure 5. Correlation between the corotation radius(i?cjj) 
scaled by semi-major axis length (ag) and the central velocity 
dispersion (V^isp) of the galaxies studied by Corsini et al (2004) 
with the TW-method. 

a clear tendency of larger velocity dispersions being associ- 
ated with faster pattern speeds (or equivalently, with smaller 
corotation radii). Large velocity dispersions introduce ex- 
tra velocity components in the velocity measurements in- 
tegrated along the line-of-sight. Since in the TW-method 
the line of sight velocity is supposed to belong to the galac- 
tic plane, the extra component introduces a systematic er- 
ror source. The systematic errors associated with th e TW 
method are discussed in detail bv lMeidt et alj ((2008). 



8 CONCLUSIONS 

In this work we showed the existence of a clear correlation 
between breaks in the metallicity distribution and corotation 
radii of spiral galaxies. To achieve this task we compiled all 
the references that we could find on rotation curves, metal- 
licity distributions, and spiral pattern speeds or corotation 
radii of galaxies. The aim was to construct a sample of galax- 
ies for which there is at least one measurement in the lit- 
erature for each one of these quantities. The references are 
presented in tables in the Appendix. Since the data from the 
literature spans many years, different authors, and different 
methods, we established criteria to standardize the data, in 
order to avoid combining data that were obtained based on 
discrepant rotation curves and/or galactic distances. 

To locate minima, inflexions or breaks in the metallic- 
ity distribution, we developed an automatic procedure which 
consists in fitting, in sequence, several low order polynomi- 
als, a pair of straight connected lines forming an angle, and 
a pair of straight lines with a step at the connection between 
them. The significance of each of these fits was evaluated us- 
ing the chi-squared by degree of freedom. Coincident values 
for minima, inflexions and breaks were registered. Almost 
all galaxies that were better fitted by other models than the 



linear fit presented coincident minima for different curves 
(different polynomial order or the set of two lines), which 
means that the minimum is not an artefact of a given fit. 

The distribution of spiral pattern speeds presents a 
small concentration around 20 km/s/kpc, so that the spi- 
ral pattern speed of our Galaxy can be considered typical 
in the context of our sample. The distribution of corotation 
radii, on its turn, presents a cutoff around 15 kpc, which is 
consistent with the spiral density wave theory, since typical 
large spiral galaxies have a radius of about 30 kpc. The outer 
Lindblad resonance is supposed to be about the boundary 
of the visible disks and corotation is expected to be about 
midway between the Lindblad resonances. We noted a minor 
systematic deviation (about 10%) of the correlation between 
metallicity breaks and corotation from the one-to-one rela- 
tion, in the sense that the metallicity breaks are at a larger 
radius. This could be due to a real shift of the break with 
respect to corotation, or to a slight error in the determina- 
tion of the corotation radius, possibly associated with the 
TW-method. 

The existence of a minimum of metallicity at corotation, 
or equivalently, of a raising gradient in external regions of a 
galaxy, as we found in many spiral galaxies presented here, 
can be explained by a star-formation recipe like SFR tx /xi? | 
Q — fi p | or similar ones. We argued in the previous section 
that this is a quite reasonable hypothesis. The presence of 
steps in the metallicity distribution is a little more difficult 
to elucidate. Our hypothesis is that it is a consequence of the 
radial flows of gas going in opposite directions on each side 
of corotation, which tend to produce a depletion of gas at 
this resonance. This creates like a barrier, with no transfer 
of gas from one side to the other, so that the two sides 
evolve independently. We do not discard the effect of other 
resonances. 

An important consequence of the existence of both the 
breaks in metallicity gradients and of metallicity steps is 
that they points towards a long term stability of the corota- 
tion radius of the host galaxies. A varying corotation radius 
would smooth out such features. In other words the spiral 
structures of many galaxies are long-lived, as discussed by 
L+8 in the case of our Galaxy, for which a minimum age of 
the spiral structure was estimated to be 3 billion years. Our 
result also favours the idea that each galaxy has a single, 
or at least a quite dominant pattern speed, otherwise the 
correlation that we found would be destroyed. 
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APPENDIX A: REFERENCES IN THE 
LITERATURE 

In Table [ATI we compile the references on corotation radii, 
spiral pattern speeds, metallicity distributions and rotation 
curves used to obtain the results in Tables [1] [2] and IA2I 
The details on the references consulted about the corotation 
radii, with standardized values for each reference can be 
found in Table IMl 
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Table A2. Reference indexes for the corotation radii (Rcr) of each galaxy in our sample. Each index is preceded by the values found 
in the papers and re-scaled according the procedures in Section 3 and the Table [T] The symbol "V" indicates two possible values and 
the symbol "-" separate a range of values. 



GALAXY 


R-Cfl [kpc ] 


IC0342 


5.8 [34]; 14.6 [125] 


NGC0224 


14.8 [125]; 16.4 [96] 


NGC0598 


3.4 [96]; 1.2 V 7.0 [82]; 5.9 [125]; 2.9 [34]; 4.3 [92]; 4.0 [21]; 4.7 [130] 


NGC0628 


4.8 [19]; 4.6 [34]; 4.6 [125]; 2.2 V 2.9 [90] 


NGC1068 


3.1 V 3.3 [94]; 1.4 [53] 


NGC1232 


12.6 [34] 


NGC1365 


13.9 [68]; 12.3 [57]; 14.1 [35]; 10 V 20.7 [124] 


NGC1530 


13.4 [124] 


NGC1566 


8.8 [23] 


NGC2403 


5.1 [125] 


NGC2543 


4.4 V 7.5 [124] 


NGC2903 


4.2 [108]; 4.0 [125] 


NGC3031 


13.1 [129]; 9.4 V 8.1 [34]; 11.9 V 12.9 [52]; 9.1 V 9.8 [70]; 11 [96]; 11.6 [97]; 12.4 [102]; 9.9 V 15.5 [110]; 

O.D [llO], y.O [lol], 1Z.O [OUJ 




A n 1 A TqoI. O 7 Ti Ql 

4.y - (A [o2J; a. 7 [loj 




E i roll a. a Tool, c 1 n oci. c o T^Ql 

o.l [oj; b.o [y»J; o.l [loo]; o.z [7oJ 


NGC3621 


6.2 [116] 


NGC4254 


13.6 [34]; 12.1 [64]; 16.5 [51] 


NGC4321 


8.9 - 7.6 V 6.4-8.9 [106]; 7.3 [20]; 8.8 [33]; 6.3 - 7.2 [95]; 7.6 V 6.8 V 8.4 [14]; 8.6 V 14 [43]; 

*7 n m n h cl. n i Tq^I. 10 roil, c o h cl. ^ 1 d c fi ncl. 7 T1 nol. n w m n To A l . 1 *? \ / 11 *? [a Al • 
/.9 - 1U.U [loj ; y.z [34J; 12 [84J; o.z [loj; o.l - 4.5 [lUbJ; ( [lUoJ; y.o V lU.y [94J; 1.3 V 11.3 [44]; 

3.3 [63]; 4.2 V 8.6 V 13.1 [55] 


NGC4736 


3.1 [50]; 0.6 [132]; 0.7 V 0.8 [78] 


NGC4826 


1.1 V 1.4 [94] 


NGC5033 


3.2 [77] 


NGC5055 


12.8 [34]; 14 [78]; 6.7 [77]; 15.4 [96] 


NGC5194 


5.7V 10.1 [33]; 5.4 V 10.2 [96]; 2.8 [136]; 5.34 V 6.7 [119]; 6.9 [46]; 2.8 - 5.6 [103]; 5.8 [34]; 
10.1 [125]; 13.2 [84]; 10.2 [110]; 5.4 [116] 


NGC5236 


3.8 [136]; 3.4 [69]; 3.4 [56]; 8.1 [115]; 4.1 [125] 


NGC5457 


1.2 [108]; 1.7 [77]; 16.7 [96] 


NGC6946 


2.8 [34]; 4.6 [125]; 1.9 [108]; 3.7 [136]; 4.9 V 9.7 [118]; 


NGC7479 


3.0 V 4.7 [123]; 10.7 [2]; 9.5 [4]; 9.5 [29]; 12.5 V 17.8 [21]; 

7.6 [93]; 14.4 [25]; 9.6 [105]; 9.3 [65]; 9.3 [91]; 0.5 [40]; 14.6 [37]; 10.0 [13] 



